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(Dated: July 19, 2005) 

Wo present a study of the decay J/ip 4). We obtain the CP-odd fraction in the final state 

at time zero, R±_ = 0.16 ± 0.10 (stat) ± 0.02 (syst), the average hfetirne of the (Bg, _b") system, 
t{B^) = 1.39lg;Jg (stat) lo'o2 (syst) ps, and the relative width difference between the heavy and 
light mass eigenstates, AF/F = {Tl - Th)/T = 0.24lo;3^ (stat) l„:o4 (syst). With the additional 
constraint from the world average of the lifetime measurements using semileptonic decays, we 
find r(B°) = 1.39 ± 0.06 ps and AF/r = Q.25to^^- For the ratio of the B° and B° lifetimes we 

obtain = 0.91 ±0.09 (stat) ± 0.003 (syst). 



PACS numbers: 13.25.Hw, 11.30.Er 



Within the framework of the standard model (SM), the mesons are expected to mix in such a way that 
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mass and decay width differences between the heavy and 
light eigenstates, AM = Mh -Ml a.nd AT = Tl - Th, 
are sizeable. The mixing phase 5(j) is small and to a good 
approximation the two mass eigenstates correspond to 
CP eigenstates. New phenomena may alter 5(j), leading 
to a reduction of the observed AF/F compared to the 
SM prediction Q. 

The decay — J/ijj (t>j proceeding through the quark 
process b — *■ ccs, gives rise to both CP-even and CP-odd 
final states. It is possible to separate the two CP com- 
ponents of the decay ^ J/ip 0, and thus to measure 
the lifetime difference, through a simultaneous study of 
the time evolution and angular distributions of the decay 
products of the J /ip and (p mesons. Angular distribu- 
tion of B° J/ip{-^ M+M") K+K-) involve three 
angles. Current statistics are such that the use of all 
three angles characterizing the final state is not yet ben- 
eficial. We use a variable particularly sensitive to sepa- 
rating the CP states, the cosine of the transversity polar 
angle (called "tranversity"), as defined below. 

The analysis of data collected with the D0 detector Q 
at the Fermilab Tevatron collider presented in this Letter 
is an extension of a recently published study Q done un- 
der the single B^ lifetime hypothesis. We perform an un- 
binned maximum likelihood fit to the data, including the 
B^ candidate mass, lifetime, and transversity, in the de- 
cay sequence B^ — > J/^ (j), J/ip A^^M > 4' ^ K^K^ . 
We extract three parameters characterizing the 5° sys- 
tem and its decay: t{B°) = l/F, where F = {Th+Tl)/2; 
AF/F; and R±, the relative rate of the decay to the CP- 
odd states at time zero. The average lifetimes of B^ and 
B°, as defined above, are expected to be equal to within 
1% and their ratio is determined by measuring the 
lifetime of in the similar decay topology of J/ tpK* . 

The data were collected between June 2002 and Au- 
gust 2004. The sample is selected by requiring two recon- 
structed muons with a transverse momentum pT > 1.5 
GeV. Each muon is required to be detected as a track 
segment in at least one layer of the muon system and to 
be matched to a central track. One muon is required to 
have segments both inside and outside the toroid. We 
require the events to satisfy a muon trigger that does not 
include any cuts on the impact parameter. The sample 
corresponds to an integrated luminosity of approximately 
450 pb-i. 

To select the B^ candidate sample, we apply the fol- 
lowing kinematic and quality cuts. Minimum values of 
momenta in the transverse plane for i?^, and K mesons 
are set at 6.0 GeV, 1.5 GeV, and 0.7 GeV, respectively. 
J/ ^ candidates are accepted if the invariant mass is in 
the range 2.90 - 3.25 GeV. Successful candidates are con- 
strained to the average reconstructed J/ip mass of 3.072 
GeV. Decay products of the (p candidates are required to 
satisfy a fit to a common vertex and to have an invari- 
ant mass in the range 1.01 - 1.03 GeV. We require the 
{J/ip,4>) pair to be consistent with coming from a corn- 
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FIG. 1: The invariant mass distribution of the (J/i/;,0) sys- 
tem for candidates. The curves are projections of the 
maximum likelihood fit (see text). 



mon vertex, and to have an invariant mass in the range 
5.0 - 5.8 GeV. In case of multiple (j) meson candidates, 
we select the one with the highest transverse momentum. 
Monte Carlo (MC) studies show that the pr spectrum of 
the (f> mesons coming from decay is harder than the 
spectrum of a pair of random tracks from the underlying 
event. We define the signed decay length of a B^ meson 
L^y as the vector pointing from the primary vertex to the 
decay vertex projected on the transverse momentum. 
To reconstruct the primary vertex, we select tracks with 
Pt > 0.3 GeV that are not used as decay products of the 
Bg candidate and apply a constraint to the average beam 
spot position. The proper decay length, ct, is defined by 
the relation ct = L^y ■ Mgo /px where M^o is the world 
average mass of the B'^ meson |^. The distribution of 
the proper decay length uncertainty a{ct) of B° mesons 
peaks around 25 /^m. We accept events with cr(ci) < 60 
/im. There are 9699 events satisfying the above cuts. 

The resulting invariant mass distribution of the 
{J/ip,4>) system is shown in Fig. ^ The curves are pro- 
jections of the maximum likelihood fit, described below. 
The fit assigns 513±33 events to 5° decay. 

Using the same data sample and analogous kinematic 
and quality cuts, we select a sample of 1913 events 
of the decay sequence B^ J/tpK*, J/tp — > /^^M^, 
K* K^TT^ , and update the B^ lifetime measurement 
reported in Ref. with larger statistics. 

We perform a simultaneous unbinned maximum like- 
lihood fit to the proper decay length, transversity, and 
mass. The likelihood function C is given by: 



N 



^ — nt/sig-^sig + (1 /sig)-^bck]i 



(1) 



where A'' is the total number of events, T];^^ (-^bck) 
the product of the mass, proper decay length, and the 
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transversity probability density functions for the signal 
(background) , and /gig is the fraction of signal in the sam- 
ple. Background is divided into two categories, based on 
their origin and lifetime characteristics. "Prompt" back- 
ground is due to directly produced J/-0 mesons accompa- 
nied by random tracks arising from hadronization. This 
background is distinguished from "non-prompt" back- 
ground, where the J /ip meson is a product of a B-hadron 
decay while the tracks forming the candidate emanate 
from a multibody decay of the same B hadron or from 
the underlying event. We allow for independent parame- 
ters for the two background components in mass, lifetime, 
and transversity. There are nineteen free parameters in 
the fit. 

For the signal mass distribution, we use a sum of two 
Gaussian functions with a fixed ratio of widths and nor- 
malizations, obtained in a fit to the signal-dominated 



subset satisfying ct/a{ct) > 5. We allow for two free 
parameters, the common mean value and the width of 
the narrow component. The lifetime distribution of the 
signal is parametrized by an exponential convoluted with 
a Gaussian function with the width taken from the event- 
by-event estimate of cr(ci). To allow for the possibility 
of the lifetime uncertainty to be systematically underes- 
timated, we introduce a free scale factor. 

The transversity distribution of the signal is deter- 
mined in the following way. The time-dependent three- 
angle distribution for the decay of untagged B'^ mesons, 

i.e., summed over and S^, expressed in terms of the 
linear polarization amplitudes and their relative 

phases 5i is 0: 



dcos9 dip dcosip 



(X 2|Ao(0)re~^ ^' cos^ ~ sin^ 6lcos^ tp) + sin^ V'd^lKO)! - sin^ 6lsin^ (p) 



+ \A^{0)\^e-' sin^ 0} + ^ sin 2^^1^0(0)117411 (0)| 005(62 - (5i)e"^ ^* sin^ sin 2ip 

v2 



V2 



|^o(0)||Ai(0)| cos52 sin 2^^ sin 26* cos ((5 - |^|| (0)||Ax(0)| cos^i sin^ V sin 261 sin 



Scf>. (2) 



In the coordinate system of the J/ -0 rest frame (where the 
(p meson moves in the x direction, the z axis is perpendic- 
ular to the decay plane of ^ K^K^ , and py{K^) > 0), 
the transversity polar and azimuthal angles {0, ip) de- 
scribe the direction of the and '0 is the angle between 
p{K^) and ^p{J /ijj) in the rest frame. 

Wc model the acceptance in the three angles by poly- 
nomials, with parameters determined using Monte Carlo 
simulations. We have used the SW.HELAMP model 
in the EvtGen generator 0, interfaced to the Pythia 
program j^. Reconstructed events were reweighted to 
match the kinematic distributions with the data. 

To obtain the one-angle (transversity) distribution, we 
integrate the three-angle distribution over the angles 
and ip. The resulting distribution depends on one free 
parameter, R±^ — |Ax(0)p. There is a small correction 
term due to the nonuniformity of the acceptance in the 
angle tp, which is proportional to |^o(0)P — J^^n (0)p. We 
use the CDF Collaboration measurement 9] of this dif- 
ference, 0.355±0.066. 

The lifetime shape of the background is described as 
a sum of a prompt component, simulated as a Gaus- 
sian function centered at zero, and a non-prompt com- 



ponent, simulated as a superposition of one exponen- 
tial for the negative ct region and two exponentials for 
the positive ct region, with free slopes and normaliza- 
tion. The mass distributions of the backgrounds are 
parametrized by first-order polynomials. The transver- 
sity distributions of backgrounds are parametrized as 
(1 + a2 cos^ + a4 cos^ 0). 

Results of the fit are presented in Figs. ^ - 01 The 
proper decay length distribution, and the transversity 
distribution, both with the fit results overlaid are shown 
in Figs. 121 and 13 Figurc^shows the one standard devia- 
tion ( one-ff) contour for cT{Bf) versus AF/F. It provides 
the best display of the uncertainty range for these cor- 
related parameters. Our best fit returns R±— 0.16±0.10 
and AF/F = 0.24t°;^8 at t{B°) = l.^9to\ll ps. 

We do a series of alternative fits, at discrete values 
of f(i?°). The results for AF/F, its one standard de- 
viation range, and the corresponding value of the like- 
lihood, are listed in Table ^ We verify the procedure 
by performing fits on a sample of approximately 50,000 
MC events passed through the full chain of detector sim- 
ulation, event reconstruction, and maximum likelihood 
fitting. We see no bias in the event reconstruction or 
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FIG. 2: The proper decay length ct of the candidates in 
the signal mass region. The curves show the signal contribu- 
tion, dashed (red); the background, lower solid line (green); 
and total, upper solid line (black) in the signal mass region. 
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FIG. 3: The distribution of the cosine of the transversity 
angle in the signal mass region, for "non-prompt" events, with 
the results of the maximum likelihood fit overlaid. 



in the fitting procedure. The fits reproduce the inputs 
(cr = 439 ^m, AT/F — 0, and a range of R± between 
and 1) correctly within the statistical precision of 2 /im 
for cr, 0.01 for Rj_, and 0.025 for AT/T. We test the 
sensitivity of the results to the parametrization of the 
signal and background mass distributions by varying the 
parameters of the two-Gaussian function. To test the 
sensitivity of the results to the background model, we 
add a quadratic term in the background mass distribu- 
tion. We find a non-negligible effect from the extra term 
in the non-prompt background on cr and AT/T. We have 
also tested the sensitivity of the results to the assump- 
tion that the lifetime and the transversity distributions 
of background are independent of mass. The effect of 
the uncertainty of the detector alignment on the lifetime 
measurement was estimated in Ref. The effects of 
systematic uncertainties are listed in Table HTl 
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FIG. 4; The D0 defauh one-cr (stat.) contour {5\n{L) = 0.5) 
compared to a one-a band for the world average (WA) mea- 
surement based on flavor-specific decays, t/s = 1.442 ± 0.066 
ps. A simultaneous fit to our data and the WA gives a one-o- 
range r(B°) = 1.39 ±0.06 ps and AT /T = 0.25to\t- The SM 
theoretical prediction is shown as the horizontal band. 



TABLE I: Fit resuhs for AF/F at fixed values of r(B^). For 
each assumed value of T{Bg), the likelihood as function of 
AF/F is symmetric and parabolic. 



r{B"s) (ps) 


AF/F 


Aln(£) 


1.23 


-0.13 ±0.15 


0.51 


1.27 


-0.03 ±0.17 


0.32 


1.31 


0.07 ±0.19 


0.17 


1.35 


0.16 ±0.21 


0.04 


1.39 


0.24 ±0.20 


0.0 


1.43 


0.31 ±0.19 


0.06 


1.47 


0.37 ±0.18 


0.20 


1.51 


0.43 ±0.18 


0.42 


1.55 


0.48 ±0.18 


0.69 



We also conduct a test with an ensemble of 1000 
pseudo-experiments with similar statistical sensitivity, 
chosen from the distribution described by Eq. ^ with 
the same parameters as obtained in this analysis. Both 
the spread of uncertainties and of the central values of 
the fit parameters are in good agreement with the results 
reported here. Our results are consistent with the CDF 
Collaboration results 0|, also shown in Fig. 01 

i?^ lifetime measurements from semileptonic (flavor- 
specific) data provide an independent constraint on the 
average lifetime and lifetime difference in the sys- 
tem. The world average ^ i?" lifetime is Tfs = 
l/F/s = 1.442 ± 0.066 ps. This resuh is based on single- 
exponential fits in the flavor-specific decay channels, 
which determine the following relation (shown in Fig. 

mi of F and AT/T: F/, = F - (AF) V2r + 0(A^)V^^ 
Applying the above constraint to our measurement, we 
obtain r(B^) = 1.39 ± 0.06 ps and AT/T = 0.25^\i 
This result is consistent with the SM expectation [ll| of 
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TABLE II: Sources of systematic uncertainty. The numbers 
reflect the variation of the fitted central values associated with 
the one-cr variation of the corresponding external input pa- 
rameters. The second item includes contributions from the 
variation of the acceptance as a function of (/? and tp, as well 
as from a one-cr variation of the quantity |j4o(0)p — |^||(0)p. 



time, we obtain 

7(B°) _ 
t(BO) ^ 



0.91 ± 0.09 (stat) ± 0.003 (syst) 



Source 


cr{B"), fj.m 


Ar/r 


R± 


Acceptance vs. cosO 


±0.6 


±0.001 


±0.005 


Integration over ip, ip 


±0.2 


±0.001 


±0.02 


Procedure test 


±2.0 


±0.025 


±0.01 


Momentum scale 


-3.0 






Signal mass model 


±1.0 


±0.009,-0.017 


±0.007 


Background mass 


-3.5 


±0.02 


-0.002 


Detector alignment 


±2.0 






Background model 


±0.5 


±0.016 


±0.005 


Total 


-5.6,±3.1 


-0.04, ±0.03 


±0.02 



0.12 ± 0.05. 

All the results presented above are obtained under a 
tacit assumption that the CP-violating phase is negli- 
gible, as predicted by the SM {Scfi = 4>ckm — —0.03). 
Future improvements on the measurement of AF/F may 
exclude models predicting large deviations of Scf) from the 
SM value. 

In summary, we have measured the CP-odd fraction for 
the decay 5° J/ip 4>, and the correlated parameters of 
the average lifetime of the (B^, B^) system r(i?°) = 1/F, 
and the relative width difference AF/F, or, equivalently, 
the mean lifetimes of the light and heavy i?^ eigenstates, 
respectively. We obtain 

i?_L =_0.16 ± 0.10 (stat) ± 0.02 (syst), 
AF/F = 0.24±0-J« (stat) 1^04 (syst), 



r(B°) = 1.39t'o,\l (stat) t^ol (syst) ps, 
TL = l.2At°o.\t (stat) (syst) ps, 

TH = 1.58t°:l (stat) t^om (syst) ps. 
We have updated the measurement of the mean life- 
time of the B'^ meson with doubled statistics. With the 
systematic uncertainty estimated in Ref. 3] , the updated 
measurement is 

t{B") = 1.530 ± 0.043 (stat) ± 0.023 (syst) ps. 
For the ratio of the average B° lifetime to the life- 
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